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We propose an order parameter to characterize valence-bond-solid (VBS) states in quan- 
tum spin chains, given by the ground-state expectation value of a unitary operator appearing 
in the Lieb-Schultz-Mattis argument. We show that the order parameter changes the sign 
according to the configuration of the valence bonds. This allows us to determine the phase 
transition point in between different VBS states accurately. We demonstrate this theory in 
the bond-alternating Heisenberg chain and in the frustrated spin ladder. 



§1. Introduction 



To study Haldane's conjecture EP, Affleck, Kennedy, Lieb and Tasaki proposed 
the valence-bond-solid (VBS) state, and showed that the ground state of the spin 
5 = 1 Heisenberg chain is described approximately by the VBS statecP. Since the 
spin configurations of the VBS state show the hidden antiferromagnetic order, den 
Nijs and Rommelse proposed the string order parameter (SOP) to characterize the 
S = 1 Haldane phase cP, 
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where a = x,y,z and |<?b) means the ground state. Thus this order parameter 
enables us to detect the VBS state indirectly. The SOP was generalized to S > 1 
cases by OshikawacP. 

On the other hand, Affleck and Lieb studied the Haldane's conjecture by the 
Lieb-Schultz-Mattis (LSM) type arguments). However, relation between the LSM 
argument and the VBS picture including the SOP has not been fully understood. In 
this paper, we show that the ground-state expectation value of a unitary operator 
appearing in the LSM argument, for a spin chain with length L, 
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plays a role of an order parameter which detects VBS ground states directly, and 
that it can be applied to determination of phase boundaries in between different VBS 
states. We demonstrate this idea in the bond-alternating Heisenberg spin chains, and 
in the S = 1/2 two-leg frustrated ladder. 



typeset using ■pTpTfX.sty <ver.l.0> 



2 



M. Nakamura and S. Todo 



§2. Properties of the Order Parameter 

2.1. The valence-bond-solid picture 

According to the LSM argument, if the unique ground state |$b) and U\&o) 
with U = ex.p[(2wi/ L) ^2 j =1 jSj] are orthogonal (zl = 0), there exists at least one 
eigenstate with energy of Now we consider the overlap zl in the opposite 

situations such as the Haldane gap phase. In order to calculate zl explicitly in VBS 
states, we introduce the Schwinger boson representation for the spin operators: 

S+ = a]bj, S-=b] aj , S^^aj-btbj), S j = ±(a]a j + b]b j ), (2-1) 

where these bosons satisfy the commutation relation [aj,a|] = [&i,bj] = b~ij with 

all the other commutations vanishing, a J (bj) increases the number of up (down) 
5=1/2 variables under symmetrization. Then, a generalized VBS state in a periodic 
system is written as 

l<B's n) ) - ^ if (4-i, 2fc ) m (< 2fc+ i) n |vac), (2-2) 

where B] ■ = a]b\-b]a\, J\f is a normalization factor, and Ivac) is the vacuum with 
respect to bosons, m and n are integers satisfying m + n = 2S. Using relations 

Uajw 1 = a]e + ' mj / L and Ub]lJ- 1 = b]e~ inj / L , a twisted valence bond UB^^U' 1 
for 1 < j < L — 1, and that located at the boundary are calculated as follows, 

UBj^U- 1 = e-^ L a]b] +1 - e^ L b]a] +1 , (2-3) 

UB^U' 1 = -{e-^ L a\b\ - e^ L b{a\). (2-4) 

In the latter case, a negative sign appears for each valence bond. Thus the asymptotic 
form of zl is given by 

zl = (^Vl^) = (-!)"[! " 0{l/L)]. (2-5) 

It turns out that zl changes its sign according to the number of valence bonds at 
the boundary, zl can be calculated in more detail by the matrix product formalism. 

2.2. The sine- Gordon theory 

Next, we consider zl from the viewpoint of a low energy effective theory. In 
general, the Lagrangian density of a spin chain is given by the sine-Gordon model, 

C = 2^K t V ^' r )] 2 - ^2 cos[V2^,r)], (2-6) 

where r is the imaginary time, a is a short range cut off, and K and y^ are the 
parameters determined phenomenologically. In the gapped (gapless) region one has 
2/$(0 — * i 00 (y</>(0 —> 0) f° r ' —> 00 under renormalization a — ► e'a. On the unstable 
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Gaussian fixed line [y^(0) = with K(0) < 4], a second-order "Gaussian transition" 
takes place between the two gapped states. In this formalism, the spin wave excita- 
tion created by U corresponds to the vertex operator exp(iv20), so that zl is related 
to the ground-state expectation value of the nonlinear term as zl oc (cos(\/2<^)) and 
the three fixed points = ±oo, correspond to z^ = =f1, 0, respectively. Thus the 
Gaussian critical point can be identified by observing zi = 

§3. Application to Physical Systems 

3.1. The bond- alternating Heisenberg chain 
In the bond-alternating Heisenberg chain, 

L 

n = J2^-s(-ms j -s j+ i, (3-i) 

the VBS picture is considered to be realized approximately: The configuration of the 
valence bonds (m, n) changes from (0, 2S) to (2S, 0) successively as S is increased 
form —1 to 1, meaning the existence of 2S quantum phase transitions. We calculate 
zl for 25 = 1,2,3,4 in a periodic system by the quantum Monte Carlo (QMC) 
method with the continuous-time loop algorithm CP , and confirmed that qualitative 
behavior of zl agrees with our theory given in Sec. 

Next, we determine the successive dimerization transition points by observing 
zl = with L up to 320. Extrapolation of the data has been done by the function 
S C (L) = 5 c (oo) + A/L 2 + B/L 4 + C/L 6 . For S = 1 [(1, l)-(2, 0)], S = 3/2, [(2, 1)- 
(3,0)] and 5 = 2 [(2, 2)-(3, 1), (3, l)-(4, 0)] cases, we obtain the transition points 
<5 C = 0.25997(3), 0.43131(7), 0.1866(7), 0.5500(1), respectively, where () denotes 2a. 
The results are consistent with the previous estimates, but much more accurate: 
5 C = 0.2595(5) for the 5=1 case by the QMC calculation for the susceptibjjitya-' , 
and 5 C = 0.2598,0.4315,0.1830,0.5505 obtained by the level-crossing method!) . 

3.2. The spin-1/2 frustrated ladder 

We discuss Zl in the 5 = 1/2 two-leg ladder with frustration: 
L 

7i> = ^[£7 • Sj+x + Tj ■ Tj + i + J±Sj ■ Tj + J x {Sj ■ Tj + i + Tj ■ Sj+i)], (3-2) 
i=i 

where Sj and Tj denote 5 = 1/2 operators at j-th site. In has been pointed 
outE§> that there appears competition between a rung dimer phase (J x <C J±) and 
a Haldane phase (Jj_ <C Jx < 1)> and that these two phases are identified by two 
distinct SOPs, Oodd'^even given by Eq.(l_l) where 5" is replaced by the following 



two types of composite spin operators: 

°odd,j = °j ' 1 j i °even,jr' = D j 1 K°'°) 

Then, dd = and O even / for the rung dimer phase, and dd an d O evcn = 
for the Haldane phase. Instead of C dd and O even , we introduce z Q dd and z even defined 
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by Eq.( |T2[ ) where S? is replaced by 

^oddj and Seven,j- In Fi S- h we snow 
z D dd arid Zeven for J± = 1 and L = 14 
calculated by the exact diagonaliza- 
tion. We find that | ^ dd I l^evenl i n 
general, and z dd and z even have oppo- 
site signs. Especially, z 0( id = z even = 
at the same point which corresponds 
to the dimer-Haldane phase bound- 
ary. This tells us that one of the 
two zl has enough information to de- 
scribe this system, since zl character- 
izes the system by its magnitude and 
sign, while the SOPs give zero or fi- 
nite value without changing the sign. 
Thus zl turns out to be a more ratio- 
nal order parameter to describe VBS 
states in a unified way. The properties 
phase boundary from numerical data c 
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Fig. 1. J x -dependence of z dd and z C ven in the 
S — 1/2 frustrated ladder for Jj_ = 1 and L = 
14, calculated by the exact diagonalization. 
The point z 0( jd = z avcn = gives the dimer- 
Haldane phase boundary. 

f zl is also useful to determine the accurate 
finite-size clusters. 



§4. Summary 

We have introduced zl given by Eq. (1-2) as an order parameter to character- 



ize VBS states. The proposed order parameter changes its sign according to the 
configuration of valence bonds. This property enables us to determine the critical 
point between different VBS states by observing zl = 0. We have demonstrated this 
theory for phase transitions in the bond-alternating Heisenberg chain, and in ±he 
frustrated spin ladder. Details of the present results will be published elsewhereliiP. 
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